Templated self-assembly of patchy particles 
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We explore the use of templated self-assembly to facilitate the formation of complex target structures made from patchy particles. 
First, we consider the templating of high-symmetry shell structures around a spherical core particle. We find that nucleation 
around the core particle can inhibit aggregate formation, a process which often hinders self-assembly. 

In particular, this new assembly pathway allows dodecahedral shells to form readily, whereas these structures never form in the 
absence of the template. Secondly, we consider the self-assembly of multi-shell structures, where the central icosahedral core is 
known to form readily on its own, and which could then template the growth of further layers. We are able to find conditions 
under which two- and three-shell structures successfully assemble, illustrating the power of the templating approach. 
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1 Introduction 

It has been a long-held goal to use self-assembly to create com- 
plex, ordered structures on the micro- and nanoscale. 1 Biolog- 
ical systems with their dazzling array of ordered and precise 
self-assembled structures 2 show what can be possible given 
sufficient control of the design and the interactions between the 
building blocks. Thus, with the goal of taking the first steps 
to achieving similar control in synthetic systems, researchers 
working on nanoparticles and colloids are making great efforts 
to synthesize "patchy" particles that have interactions only in 
specific directions. 2«ES 

The formation of icosahedral virus capsids, proteinaceous 
shells with specific size and structure that are designed to en- 
capsulate the viral genome, provides one of the archetypal ex- 
amples of biological self-assembly, and also one of the most 
studied. 11 If similar structures are to be achieved in synthetic 
systems, it will be important to understand the basic physi- 
cal principles of such self-assembly and the design rules for 
the interactions between the constituent particles. So far, the- 
or yE2IIH anc [ simulations^^ have been most concerned with 
understanding the assembly of smaller (T = 1, 3 or 4) empty 
capsids. However, functional viruses are not empty, but contain 
the genomic material. For icosahedral RNA viruses, the capsid 
proteins and the RNA typically co-assemble, and experiments 
have indicated that the kinetics of such assembly can be differ- 
ent from that of empty capsids. 29 Indeed, this ability of capsid 
proteins to encapsulate is not limited to the viral genome, but 
has been exploited to achieve capsid assembly around nanopar- 
ticles,^!!! nanoemulsion droplets 33 and anionic polymers™ 
Furthermore, for larger and more complex capsids, assembly 
can be dependent on the presence of "scaffolding" proteins, 
wh ich ar e thought to template the correct assembly of the cap- 

^ [35l36l 
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Although less studied, simulations'^!!^! a nd theory^E^have 
begun to address the assembly of filled capsids. Particularly 
relevant to the present study are the simulations of Hagan and 
coworkers, who have shown that the change in assembly mech- 
anism from homogeneous nucleation of the capsid proteins in 
the empty capsid case to heterogeneous nucleation around a 
core enhances the ability of the capsids to assemble.^ 

When considering the design of synthetic systems of patchy 
particles which assemble into capsid-like structures, one has to 
consider the potential differences in the interactions between 
the biological and synthetic cases. For capsids, the interfaces 
responsible for the protein-protein interactions not only have 
to be in contact, but have to have the correct relative orienta- 
tion. However, first-generation patchy colloids and nanopar- 
ticles are unlikely to have this "torsional" component in the 
potential. Simulations have indicated that the one-component 
self-assembly of simple monodisperse targets is still feasible in 
the absence of torsional constraints;^^! however, there can be 
significant differences in the mechanisms of assembly and the 
nature of the kinetic traps compared to virus capsids. In partic- 
ular, the lack of a torsional component in the interparticle po- 
tential to enforce convexity in the growing clusters le ads to dis- 
ordered aggregation competing with correct assembly .^^Fur- 
thermore, as the size of the target structure increases, the diffi- 
culty of assembly increases much more rapidly when torsional 
constraints are not present. For example, the self-assembly of 
20-particle dodecahedra occurs readily with a protein-like po- 
tential with torsional constraints,^ but is seemingly impossible 
without. 46 

Here, we will explore whether templated assembly might po- 
tentially provide a way to allow the formation of more complex 
targets using synthetic patchy particles. We will consider both 
assembly around a single central core particle (Section [3} and 
also the assembly of two- and three-shell structures (Section 
[4]), where the central shell is known to readily assemble on its 
own J15] if me ma terials for the different particles were chosen 
appropriately, the template could then be selectively removed 
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by chemical 47 or thermal 48 treatment, if desired. 

We should also note that templating has been a useful strat- 
egy in supramolecular chemistry, 49 and even for the macro- 
scopic assembly of magnetic subunits into polyhedral shells. 

2 Methods 
2.1 Potential 

To model the patchy particles, we use the pot ential that we have 
used in our previous work on self- assembly,! 45 1 46 1 but with the 
additional feature that we consider multiple particle types. The 
model has also been used to study the crystallization of patchy 
colloids, EES a nd, with an additional torsional component to 
the potential, the self-assembly of protein complexes. 

In the model, the repulsion between particles i and j is based 
upon an isotropic Lennard- Jones potential 
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but where the attraction is modulated by an orientational term, 
V ang . Thus, the complete potential is 



where Qi is the orientation of particle i, and 
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Fig. 1 Target geometries: (a) AB12, a centred icosahedron and (b) 
AB20, a centred dodecahedron. 



2.2 Simulations 

We use the virtual move Monte Carlo algorithm of Whitelam 
and GeisslerE3SSI to simulate the dynamics. We choose this al- 
gorithm for a number of reasons. Firstly, as with other Monte 
Carlo algorithms based on local moves, the algorithm gives dif- 
fusive dynamics, as is appropriate for a model of colloids and 
nanoparticles in solution. Secondly, it can generate collective 
motion of the particles. Thirdly, the algorithm is designed to 
generate the correct relative diffusion rates for clusters of dif- 
ferent size. The latter two features represent potential advan- 
tages over the single -particl e Monte Carlo we typically used in 
our previous studies. However, as the main mode of clus- 
ter growth in those studies was by monomer addition, this was 
not a significant drawback. Here, however, we felt that it was 
important not to artificially disfavour growth by the addition 
of clusters particularly in the case of the multi-shell-structures 
studied in Section |4j As the actual algorithm is quite involved, 
we do not give the details here, but refer the interested reader to 
the original papers J2ED We note that we use a maximum trans- 
lational move size of 0.3 a m i n (where a m i n is a for the smallest 
particle in the system), as we have found this to give a good bal- 
ance between moves that lead to internal rearrangements within 
clusters and moves that lead to diffusive behaviour of complete 
clusters.^ 



where ai j is the angle between the normal to patch a on parti- 
cle i and the interparticle vector r^ , and the 'max' selects the 
pair of patches that have the strongest interaction for the current 
geometry. We assume that the particle sizes are additive, i.e. 

crij = (era + c r jj)/2. We also generally choose a pw , a measure 3 Assembly around a core particle 
of the width of a patch, to be the same for all patches, except in 
the case of the central particle in Section [3] By contrast, we al- 
low the well-depth of the patch-patch interactions, e a p to vary 
(e re f = max[^]). In addition, for computational efficiency, 
we cut and shift the potential at 3.5 cr max (where cr max is a for 
the largest particle in the system), and also shift the crossover 
distance in Eq. ^ so that it still occurs when the potential is 
zero. 

We note that in this model, we have effectively coarse- 
grained out the solvent degrees of freedom, and so when we 
talk about a gas phase in our model, this would correspond to a 
dilute solution. 



In this section, we consider the growth of symmetric shell-like 
clusters around a spherical core particle A. The two cases that 
we consider are of an icosahedral and a dodecahedral outer 
shell, as depicted in Fig. [T] The B particles that make up 
the outer shell have m identical patches whose positions are 
such that they point directly at the neighbouring B particles 
in the target geometry. These patches only interact with the 
equivalent patches on other B particles and with an interaction 
strength sbb • The B particles also have a second type of patch 
which points directly towards the centre of the target clusters. 
These patches only interact with the A particles and have an 
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interaction strength e A B- For this AB interaction, there is no 
dependence on the orientation of particle A - its contribution 
to the attraction is isotropic (or equivalently, l/a pw ^A = 0). 

Figure [2] shows the results for the assembly of centred icosa- 
hedra as a function of patch width and temperature for two dif- 
ferent values of Sab • For comparison, the yield of icosahedra in 
simulations with only B particles present is shown in Fig.[2jc). 
In this case, which we have considered in detail previously, 45 
there are two basic mechanisms of assembly. Firstly, at temper- 
atures close to the clustering temperature, T c , at which icosahe- 
dral clusters become stable with respect to a gas of monomers, 
and at smaller a pw , assembly of the icosahedra proceeds by 
direct nucleation. By contrast, for wider patches and lower 
temperatures, large aggregates form first, but these can then 
undergo further internal rearrangements leading to the forma- 
tion and budding off of complete icosahedra. However, at even 
lower temperatures, the time scale for the rearrangements of the 
aggregates becomes so long that no icosahedra are able to form, 
and at even wider patch widths liquid aggregates are thermody- 
namically stable with respect to icosahedra. 

In our simulations with both A and B particles present, there 
is the potential for competition between the formation of cen- 
tred and uncentred icosahedra. Interestingly, we find that for 
£ab/^aa = 0.5 (Fig. [2ja)), centred icosahedra are the domi- 
nant product in the region of parameter space that we previ- 
ously identified as being dominated by direct nucleation for 
the pure B system, and that uncentred icosahedra preferen- 
tially form in the region dominated by the budding-off mech- 
anism. A closer comparison shows that the centred icosahe- 
dra start to form at a somewhat higher temperature than for 
the pure B system, which is unsurprising due to the extra sta- 
bilization due to interactions with the central particle. In the 
region T c (Bi2) < T < T c (ABi2), there is no competition be- 
tween the two forms as the only cluster stable with respect to 
the monomeric gas is the centred icosahedron and so growth 
is expected to occur by templated assembly around the central 
particle. 

As one moves below T c (Bi 2 ), the uncentred icosahedra can 
also start to form, and in the region where aggregation is ini- 
tially more rapid than cluster formation, the uncentred icosa- 
hedra preferentially form. This preference arises because the 
formation of liquid aggregates is driven by the BB interactions 
and the A particles are generally excluded from the interior of 
these aggregates. Again, at too low temperatures, the system 
gets trapped in aggregates, rather than forming clusters. 

For £ab/£bb = 1, (Fig.[2jb)) the increased stabilization of 
the centred relative to the uncentred icosahedra is such that the 
centred icosahedra start to form at significantly higher temper- 
atures than for the pure B system, and the uncentred icosahedra 
rarely form in any part of the parameter space. 

The effect of £ab/^bb is explored further in Figure [5] where 
we look at the dependence of the competition between centred 
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Fig. 2 (a) and (b) Yields (averaged over five repeats) of centred and 
empty icosahedra, and the number of particles in aggregates after 10 6 
MC cycles as a function of patch width (measured in radians) and the 
temperature for a system of 120 B particles and 10 A particles at a 
density of B particles of 0.15 cr^~ B 3 , where ctaa/obb = 0.95 and (a) 
£ab/£bb = 0.5 and (b) £ab/£bb = 1.0. For comparison in (c) we 
have plotted the yield of icosahedra and the number of particles in 
aggregates under identical conditions but where only B particles are 
present. Clusters are identified as centred and empty icosahedra if 
they have the correct number of particles and within two of the 
expected number of bonds, i.e. 28-30 and 40-42 bonds for empty and 
centred icosahedra, respectively, where an interaction is considered a 
bond if it is at least 40% of the well depth. Aggregates are defined as 
clusters containing at least 25 particles. 
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Fig. 3 The yields of centred and empty icosahedra after 10 6 MC 
cycles as a function of sab/sbb and temperature for a system of 120 
B particles and 10 A particles at a density of B particles of 0.15 cr^" B 3 , 
where a pw = 0.35 and ctaa/ctbb = 0.95. Definitions of icosahedra 
and aggregates are the same as in Figure [2] 



and uncentred icosahedra on s A b / £bb at a value of the patch 
width for which the yield of centred icosahedra is high in Fig. 
[2] (a, b). At £ab/^bb = 1, there is a wide range of tempera- 
ture over which centred icosahedra successfully assemble, and 
the templated assembly of the icosahedron around A particles 
dominates over direct nucleation of Bi 2 icosahedra. Below this 
temperature window, the system forms kinetic aggregates (at 
this patch width, the local structure of these aggregate is not 
that similar to the target, so little product results from rear- 
rangement of these aggregates, unlike at larger patch widths). 

As £ab/£bb decreases, T c (ABi2) decreases, but the tempera- 
ture at which aggregation begins remains relatively unchanged, 
because this is mainly determined by Ebb- Hence, the win- 
dow over which successful AB12 assembly occurs decreases. 
Furthermore, at sufficiently small ^ab/^bb, T c (ABi2) becomes 
lower than T C (B 12) and at this point the energy gained from 
an A atom being inside the Bi 2 icosahedron does not offset 
the loss of entropy. A simple estimate of the value of Sab / £bb 
for this crossover can be found using the approximation that 
T c oc E gs / (n — 1) where E gs is the ground state energy of the 
cluster and n is the number of atoms in the cluster; this ex- 
pression has been found to provide a surprisingly accurate de- 
scription of the dep endence of transition temperatures on rel- 
ative patch strengthJ^EED Using E gs (B 12 ) « — 30<Sbb and 
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-30£bb — 12<sab gives a crossover value of 



Sab/sbb = 5/22 = 0.227. 

However, even before this value of £ab / £bb is reached, the 
mechanism of assembly of AB12 will have started to change. 
At Sab = £aa, ^(AB n _i) < E(B n ) for any n and so it 
is favourable for the B particles to grow around the templat- 
ing A particle. However, as £ab/^bb decreases, the value of 
n at which it becomes more energetically favourable to form 
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Fig. 4 The percentage yields of (a) centred icosahedra and (b) 
13-mers. 10 6 MC cycles as a function of ctaa/ctbb and temperature 
for a system of 120 B particles and 10 A particles at a density of B 
particles of 0.15 (jg B 3 , where a pw = 0.35 and £ab/£bb=1. The 
definitions of icosahedra and aggregates are the same as in Figure [2] 
and 13-mers are simply clusters with 13 particles. 



AB n _i rather than B n increases, and the role of A as a tem- 
plate diminishes. Consequently, at £ab/£bb = 0.227, the direct 
nucleation of uncentred icosahedra already dominates. 

So far, we have considered cases where the A particle is the 
right size to fit inside an icosahedron. In Figure]?] we consider 
the effect of the size of the A particle on the ease with which 
centred icosahedra form. It can be seen that there is a lim- 
ited size range over which the centred icosahedra form. When 
caa/obb is too large, although templated growth of AB n clus- 
ters will occur, the curvature of the central particle is too small 
to allow the particles to form icosahedra. When oaa/cbb is 
too small, although the initial nucleation of the icosahedra still 
occurs on the A particle, not all the particles in the growing 
icosahedra can maintain contact with the A particle. Thus al- 
though AB12 clusters still result, the A particle is now in an 
off-centre position touching only a subset of the B particles. 
Figures |4ja) and (b) differentiate between the yields of AB12 
clusters with the A particle centred or off-centred. It is no- 
ticeable that the off-centred clusters persist to lower oaa/cbb 
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but that the temperature window over which they can form de- 
creases with decreasing oaa/cbb, because the energetic sta- 
bilization of the cluster provided by the A particle decreases, 
as it can contact fewer and fewer of the B particles. Only at 
the smallest values of oaa/cbb do unfilled icosahedra begin to 
form. 

Although icosahedral clusters can readily assemble in our 
simulations with or without a templating central particle, the 
same is not true for other target clusters. In particular, we 
previously found that it was impossible to get appropriately- 
designed patchy particles to form 20-particle dodecahedral 
shells. 46 The essential problem is that the system always 
prefers to form aggregates rather than clusters, because the ag- 
gregates are first to become thermodynamically stable as the 
temperature is decreased, i.e. T aggreg > T c (£? 2 o)> and so there 
is no temperature window for which the target clusters are the 
only species stable with respect to the gas. Furthermore, even 
when the dodecahedral clusters are more stable, i.e. for T < T c 
the aggregates form more rapidly than dodecahedral clusters. 

Here we investigate whether we can get dodecahedral clus- 
ters to form using templated self-assembly, the idea being that 
the addition of templating particles could stabilize the dodeca- 
hedra sufficiently such that T c (AB 2 o) > ^kggreg* thus resolv- 
ing the thermodynamic problem noted above. Moreover, tem- 
plated growth around the central particles will help kinetically, 
by forcing the clusters to grow with the correct curvature. 

It can be seen from Figure [5ja) that there is now a clear re- 
gion of parameter space where centred dodecahedra form. As 
expected, no uncentred dodecahedra form, and so the compe- 
tition is simply between templated assembly and aggregation. 
The role of Sab in stabilizing the target structure and enabling 
assembly is clear from Figure[5] which shows the yield of AB20 
dodecahedra as a function of Sab / £bb • In the temperature win- 
dow T aggreg < T < T c (AB2o), dodecahedra now readily form, 
as the target clusters are the only species stable with respect 
to the gas. Below this window, aggregation dominates because 
it occurs more rapidly than cluster formation, and because the 
aggregates are so structurally different from the target that rear- 
rangement of the aggregates to form the target will never occur. 

As £ab/£bb decreases, there is a corresponding decrease in 
T c , and so the window of successful assembly first narrows and 
then disappears at £ab/£bb = 0.15, where T c « T aggreg . Be- 
low this value, similar to the pure B system, dodecahedra are 
never found to assemble and aggregation always dominates. 
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Fig. 5 The yield of centred and empty dodecahedra and the number 
of particle in aggregates (a) as a function of patch width and 
temperature at sab/sbb = 1.0 and (b) as a function of sab/sbb and 
temperature at cr pw = 0.35. The system consists of 10 A particles and 
200 B particles at a density of B particles of 0.15 cr^" B 3 , where 
obb/ctaa = 1-80. The simulations were of length 10 6 MC cycles 
(210 x 10 6 MC steps). Clusters are identified as centred and empty 
dodecahedra if they have the correct number of particles and within 
two of the expected number of bonds, i.e. 28-30 and 48-50 bonds for 
empty and centred dodecahedra, respectively. Aggregates are defined 
as clusters containing at least 41 particles. 



4 Multi-shell assembly 

The above results for the dodecahedra illustrate the potential 
role of templated self-assembly in enabling the assembly of 
more complex target structures in systems of patchy particles 
without torsionally- specific interactions. In this section we take 
this one step further to consider the self-assembly of multi-shell 



5 




Fig. 6 Multi-shell target structures: (a) A12B20, a dodecahedron 
surrounding an icosahedron, and (b) A12B20C30, an 
icosidodecahedron surrounding a dodecahedron which itself 
surrounds an icosahedron. There are two views of each structure, the 
ones on the right having the outer shell partially transparent to give a 
clearer view of the core. The clusters were produced by 
self-assembly at T = 0.14 £aa& -1 , and so show some deviations 
from the ideal structure due to thermal vibrations. The relative 
particle sizes were chosen to allow the structures to form without any 
strain, i.e. ctrb = 1.213 ctaa and ace = 1.851 ctaa- 



structures, where the outer shell(s) surround a polyhedral core 
that we already know to assemble readily. 

Similar highly- symmetric multi- shell structures are seen in 
biology. For example, some icosahedral viruses have capsids 
with two or more proteinaceous sh ells, either in their native 
stated or in malformed structures J 60 * 61 1 Furthermore, multi- 
enzyme complexes have been discovered that have an open 
multi-shell structure, e.g. in some species the pyruvate dehy- 
drogenase complex has a dodecahedral inner core of 20 El 
trimers which can be surrounded by 60 E2 tetramers^ or E3 
dimers^to give a complex with overall icosahedral symmetry. 
However, the particles that we consider, unlike these proteins, 
do not have the advantage of torsionally-specific interactions. 

The two example structures that we consider are illustrated 
in Fig. [6] At the centre of both is an icosahedron of 12 A par- 
ticles. In the first example, this icosahedron is surrounded by 
a dodecahedron of 20 B particles. In the second example, in 
addition to the dodecahedral shell, there is a further outer shell 
of 30 C particles with the structure of an icosidodecahedron. 

The general principle for choosing the geometry of the 
patches for these particle is that in the ideal target cluster there 
will be patches pointing directly at the neighbours both in the 



current shell and in the adjacent shell(s). So, the A particles 
will have five patches that interact with strength eaa with the 
equivalent patches on other A particles and that are responsible 
for forming the icosahedron. Each A particle will also have a 
further 5 patches that point at the pentagon of B particles that 
surrounds each vertex of the icosahedron in the target structure. 
These latter patches interact with strength Eab with the three 
patches on the B particles that point to the three A particles of 
the triangular face of the icosahedron that each B particle sits 
above in the target cluster. Finally, the B particles also have 
a further three patches that interact with strength ebb with the 
equivalent patches on other B particles and that are responsible 
for forming the dodecahedron. Thus, the ground state energy of 
the A12B20 cluster is approximately — 30£aa — SO^ab — 30£bb- 
For the A12B20C30 cluster, the B particles have an additional 
three patches that point at the triangle of C particles that sur- 
rounds each vertex of the dodecahedron in the target structure. 
These latter patches interact with strength sbc with the two 
patches on the C particles that point to the two B particles of the 
edge of the dodecahedron that each C particle sits above in the 
target cluster. The C particles also have a further four patches 
that interact with strength sqc with the equivalent patches on 
other C particles and that are responsible for forming the icosi- 
dodecahedron. Thus, the ground state energy of the A12B20C30 
cluster is approximately -30e A A - 60£ A b — 30£ B b - 60e B c - 
60£ C c. 

The number of potential parameters when considering the as- 
sembly of these multi-shell clusters is considerably larger than 
for the structures in Section [3] Therefore, we always keep the 
size of the particles fixed at their ideal values (as given in Fig. 
[6]), and for all patches a pw = 0.4. We always use a volume 
fraction of 0.0785, which is equivalent to a number density of 

0. 15 a]jl in a one-component system of A particles. We then 
consider the effects of varying the interaction strengths of the 
patches on the self-assembly behaviour. 

Results for the self-assembly Ai 2 B 2 o are considered in Fig. 
[7] as a function of eab and sbb for two temperatures. At the 
first temperature we consider, T = 0.14 £ A a& -1 < T c (Ai 2 ), 

1. e. icosahedral A12 clusters are stable irrespective of the val- 
ues of Sab and Ebb- Thus, in the bottom left-hand corner of 
Figs. |7ja)-(c) corresponding to small Eab and £bb, the forma- 
tion of isolated A12 clusters is observed. As either eab or ebb 
is increased, A12B20 clusters become stabilized and a region 
of parameter space is reached where it becomes favourable to 
form these clusters. Close to where they first become stable, 
A12B20 clusters readily form, with yields of 80% or more com- 
mon. In this region, the majority of A12 clusters are part of 
complete A12B20 clusters and provide a stable intermediate for 
the templated growth of the target cluster. The mechanism of 
assembly is likely to be hierarchical with A12 clusters forming 
first, followed by the growth of the second shell by the addition 
of B monomers. 
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Fig. 7 The dependence of the self-assembly of A12B20 on sab and sbb at a given temperature: (a)-(c) T = 0.14 £aa& -1 , and (d)-(f) 
T — 0.18 £aa& -1 . (a) and (d) show the yield of A12 icosahedra (irrespective of whether they are bonded to B atoms), (b) and (e) show the 
yield of Ai 2 B 2 o, and (c) and (f) show the average cluster size, all after 625 000 MC cycles. Each simulation contains 120 A and 200 B 
particles, so that a maximum of ten target clusters could be formed. The volume fraction is 0.0785. 
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Fig. 8 The dependence of the self-assembly of A12B20C30 on ebc and sec at (a)-(c) T — 0.12 £aa& _1 , £ab/£aa = 0.5 and sbb/^aa = 0.7, 
and (d)-(f) T — 0.18 £aa& -1 , £ab/^aa = 1.1 and sbb/^aa = 1.1. (a) and (d) show the yield of A12B20 (irrespective of whether they are 
bonded to C atoms), (b) and (e) show the yield of A12B20C30 and (c) and (f) show the average cluster size, all after 625 000 MC cycles. Each 
simulation contains 120 A, 200 B and 300 C particles, so that a maximum of ten target clusters could be formed. The volume fraction is 
0.0785. 
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However, as £ab, and to a lesser extent £bb, is increased 
further, malformed structures also become stable with respect 
to the gas phase and the yield of both A12 cores and A12B20 
clusters falls off dramatically. It is noti ceabl e that, unlike the 
one-component self-assembling systems " 5 * 46 ! and even the sys- 
tems considered in Sect. [3] this fall-off in the yield is not as- 
sociated with the formation of system- spanning aggregates, but 
instead the average cluster size remains similar to that of the 
target cluster. Thus, there is a major difference in the configu- 
rations responsible for kinetic trapping and this is because the 
trapping is dominated by interactions between different particle 
types rather than interactions between like particles. For exam- 
ple, we can see from Fig.|5]that for aggregation between just B 
particles, kT /ebb must be below approximately 0.07, which at 
T = 0.14 £ A a& -1 corresponds to £bb/^aa > 2. Kinetic trap- 
ping is instead caused by the rapid binding between the A and 
B particles, before the A particles have assembled into icosa- 
hedra, due to the strength of the AB interactions, which in turn 
interferes with assembly for a number of reasons. Firstly, there 
are kinetic effects. The B particles make it more difficult for 
the A particles to which they are bound to come together, both 
because they slow down the diffusion of the clusters and be- 
cause the steric effect of the B particles means that it is less 
likely that a collision will lead to binding. Furthermore, in col- 
lisions between clusters, these clusters would also have to have 
'matching' structures for them to be able to come together to 
generate a cluster that retains the correct structure (i.e. that is 
still a sub-cluster of the target structure). Secondly, there is the 
thermodynamic problem that, due to the stabilizing effects of 
the AB interactions, clusters with incorrectly formed structures 
for the A core are now stable with respect to the gas phase. Al- 
though the clusters formed have on average a similar size to the 
target clusters, the system is typically a mixture of larger clus- 
ters with of the order of the 60 particles and smaller clusters 
with less than 10 particles. 

At the second temperature that we consider, T = 0.18e/c _1 , 
isolated A\ 2 icosahedra are unstable, and so can only form due 
to the stabilization provided by the second shell of B parti- 
cles. Hence, in the bottom left-hand corner of Fig. [TJd) and 
(f) the system is now a mainly monomeric vapour. Similar to 
the results for the lower temperature, as eab and Ebb are in- 
creased, there is a band in Fig. [TJe) associated with reliable 
self-assembly of the Ai 2 B 2 o target, before the yield again falls 
off due to kinetic trapping. However, this band occurs at larger 
values of eab and £bb as compared to Fig.[7Jb) because stronger 
AB and BB interactions are required to compensate for the 
higher temperature. Furthermore, as A12 clusters are no longer 
a stable intermediate, the mechanism of assembly must be more 
cooperative. Indeed, it is noticeable from a comparison of Figs. 
[TJd) and (e) that virtually all A i2 clusters that form do so as 
part of the target structure. 

These results show that the templating strategy is again suc- 




Fig. 9 Kinetic aggregates formed for the A12B20C30 target when 
T = 0.12£aa/c _1 , £ab/£aa = 0.5, £bb/£aa = 0.7, £ B c/e A A = 1.2, 
and^cc/^AA = 1.6. 



cessful in leading to the formation of structures (dodecahedral 
shells) that are otherwise impossible to form. In the next ex- 
ample, A12B20C30, we take this a stage further to show that 
templating can be used to create even more complex struc- 
tures. Given the large number of parameters for this system, 
we choose Eab and Ebb from within the region where Ai 2 B 2 o 
clusters were found to reliably form at the relevant tempera- 
ture in the previous example. In Figure [5] we then show the 
self-assembly behaviour as a function of ebc and sqc at two 
different temperatures. 

The self-assembly behaviour of this system is broadly simi- 
lar to the that for the two- shell target. Again there is a diagonal 
band of successful assembly as a function of the two interaction 
strengths (Fig.[8jb) and (e)) with the interactions insufficient to 
stabilize the target at low £bc and ecc, but so strong that the 
system becomes kinetically trapped in incorrect configurations 
at high sbc and e C c- 

However, there are also a number of differences between the 
two systems. Firstly, the maximum yields are somewhat lower 
(about 70%) but this is unsurprising given the greater complex- 
ity of the target. Perhaps more surprising is that it decreases 
by so little; this is testament to the robustness of the templated 
self-assembly approach. Secondly, there is a stronger depen- 
dence of the behaviour on £qq than there was on £bb in the pre- 
vious example, but this is simply because there are four such 
patches on each C particle (rather than the three for B parti- 
cles). Consequently, Sqc plays a greater role in the stability of 
the target structure (but also malformed structures). Thirdly, at 
large £qc the formation of large aggregates is now found to oc- 
cur. Because of their four self-interactions, C particles can start 
to aggregate at higher values of kT /sec- An aggregate from 
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this region is illustrated in Fig. [9] Binding of B particles to the 
aggregates of C particles prevents the formation of any A12B20 
clusters, but A12 icosahedra can still form because they have no 
interactions with the C particles. 

There is also an interesting difference between the self- 
assembly behaviour of the systems at the two temperatures con- 
sidered. Although for both temperatures Eab and £bb were cho- 
sen with the aim of making A12B20 clusters stable, this seems 
to have been only partially successful for the higher tempera- 
ture with yields of only 30-40% in the bottom left of Fig.^d). 
Presumably, the lower yield is partly related to the fact that we 
have chosen to keep the overall volume fraction of particles the 
same in the two examples, meaning that there is a lower con- 
centration of A and B particles compared to the two-shell case, 
and hence a lower driving force for Ai 2 B 2 o formation. Conse- 
quently, at T = 0.12£aa& -1 , the assembly mechanism can be 
hierarchical with Ai 2 B 2 o clusters a stable intermediate, lead- 
ing to a relatively broad band of high yield for the target in Fig. 
[8jb). By contrast, at T = 0.18£aa& -1 a more cooperative as- 
sembly mechanism is required to achieve high yields, because 
of the lower stability of A12B20 clusters. In this case, as the 
interaction strengths increase there is initially a broad band of 
weak assembly of the target in Fig. [SJe) with only the 30-40% 
of particles that are able to form A12B20 clusters going on to 
form the target. Only at higher interaction strengths is there a 
narrow band of higher yield, presumably because the free en- 
ergy barrier to direct nucleation of the A12B20C30 target is now 
lower. 

5 Conclusions 

In this paper, we have used computer simulations to investi- 
gate the efficacy of a templating strategy to facilitate the self- 
assembly of high symmetry monodisperse shell structures from 
model patchy particles. Importantly, in contrast to the inter- 
actions between proteins, these particles lack a torsional com- 
ponent to the patch-patch interactions, thus providing a model 
for the synthetic patchy colloids and nanoparticles that many 
groups are seeking to develop, and allowing us to explore what 
structures it might be possible to assemble with such particles 
if the patch positions and the interaction strengths could be 
precisely controlled. Previous work has shown that although 
simple target structures, e.g. 12-particle icosahedra, can readily 
form, it is hard to form more complicated structures, because 
the lack of torsional specificity in the interactions means that 
the structure of the growing clusters is not tightly controlled. 
Consequently, there is typically a competition between correct 
assembly and the formation of disordered aggregates. 

The potential advantage of using a templating strategy is that 
it opens up a new assembly mechanism, namely heterogeneous 
nucleation around the template rather than direct homogeneous 
nucleation. The simple icosahedral example allowed us to ex- 



plore under what conditions the templating pathway can domi- 
nate. Templating is generally more successful under conditions 
away from where aggregation offers a competing pathway, e.g. 
higher temperatures and narrower patch widths. Furthermore, 
the region dominated by templated assembly can be enhanced 
by increasing the interaction strength between the template and 
the assembling particles, and thus opening up a larger tempera- 
ture window over which the target structure is the only species 
stable with respect to the monomeric gas. 



We have then demonstrated the potential for templating to 
aid the formation of more complex structures by assembling 
dodecahedral and multi- shell clusters. The dodecahedral ex- 
ample is particularly noteworthy, because without the template, 
dodecahedra are never able to form as aggregation always dom- 
inates over assembly. 46 However, the stabilization of the target 
structure by the template allows a temperature window to be 
opened up where the centred dodecahedra are the only species 
stable with respect to the monomeric gas, and in this region 
the dodecahedra can now assemble relatively easily. In addi- 
tion to this thermodynamic effect, in removing aggregates as a 
competing state, templating also aids the dynamics of assem- 
bly by helping clusters to grow with the correct curvature and 
structure, even though their interactions are not torsionally spe- 
cific. For example, although the angles between the patches 
allows dodecahedron-forming particles to form hexagonal as 
well as pentagonal rings, the former are likely to be disfavoured 
because they are less congruent with binding to the template. 
Similarly, if as a cluster grows around a template, it incorpo- 
rates some kind of defect, the propagation of this defect during 
further growth of the cluster is likely to inhibit the binding to 
the template and hence reduce the stability of the growing clus- 
ter, making the annealing out of that defect more likely. Never- 
theless, even given these advantages, it is impressive that tem- 
plated assembly can allow us to form such complex structures 
as the three-shell A12B20C30 cluster in relatively high yields. 



In our simulations, we are straightforwardly able to increase 
the number of particle types and types of patches, and to intro- 
duce specificity into the patch-patch interactions. By contrast, 
even though the synthetic strategies for producing patchy par- 
ticles are rapidly improving, some of the particles whose be- 
haviour we have analysed here would be very challenging to 
synthesise, particularly in terms of the control of patch position 
and identity — DNA-mediated interactions provide a potential 
route to achieve the required specificity in the patch-patch in- 
teractions.^ However, our results also indicate that the use of 
templated self-assembly provides a means to greatly increase 
the repertoire of structures into which such particles could as- 
semble. 
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